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PREFACE 

The  research  on  this  program  was  conducted  by  the  Radar  Applications 
Division,  Radar  and  Instrumentation  Laboratory,  Engineering  Experiment 
Station,  Georgia  Institute  of  Technology,  Atlanta,  Georgia,  with  Mr.  J.  L. 
Eaves  serving  as  Project  Director,  and  Mr.  Donald  J.  Lewinski  serving  as 

Principal  Investigator.  "  - - 

This  report  was  prepared  by  the  Engineering  Experiment  Station  at 
the  Georgia  Institute  of  Technology  for  the  U.S.  Air  Force,  Rome  Air  De¬ 
fense  Center,  under  Contract  F330602-78-C-012Q.  The  Technical  Monitor 
was  Mr.  William  L.  Simkins  and  Mr.  J.  T.  Mas  t  rang  el  o  was  the  Contracting 
Officer.  For  the  purposes  of  internal  control  at  Georgia  Tech,  the  ef¬ 
fort  was  designated  Project  A-2199.  The  final  report  summarizes  the  work 
perform®*!  under  thi3  contract  and  gives  recommendation  for  future  work. 
This  report  covers  work  which  was  perforated  from  15  June  1978  to  30  Sep¬ 


tember  1978. 
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I.  INTRODUCTION 


Radar  analysis  and  system  design  are  usually  performed  under  the  assump¬ 
tions  that  radar  clutter  is  stationary  and  homogeneous;  i.e.,  the  probability 
density  and  moments  (statistics)  of  the  clutter  power  are  assumed  to  be  con¬ 
stant  in  time  for  a  single  resolution  cell  (stationary)  and  constant  in  space 
from  resolution  cell  to  resolution  cell  (homogeneous) .  In  practice  there  are 
many  situations  which  do  not  satisfy  these  conditions.  For  example,  a  single 
clutter  cell  composed  of  a  single  type  of  vegetation  which  is  moving  in  the 
wind  will  have  a  first  order  density  function  and  moments  which  vary  with 
time  and  wind  velocity.  A  radar  viewing  the  ocean  with  sufficient  resolution 
to  resolve  a  single  wave  will  record  time  varying  density  functions  and  mean 
values  for  a  single  resolution  cell,  depending  on  whether  a  wave  is  in  the 
cell  or  not.  A  single  cell  composed  of  several  different  features  (grass, 
trees,  water,  etc.)  will  have  different  statistics  than  a  cell  composed  en¬ 
tirely  of  a  single  feature.  Even  when  a  resolution  cell  is  composed  of  a 
single  feature  such  as  grass  or  trees,  the  variation  in  moisture  content, 
number  of  scatterers  per  unit  area,  etc.  result  in  an  uncertainty  in  the  den¬ 
sity  function  and  statistics.  The  feature  or  mixture  of  features  in  a  resolu¬ 
tion  cell  may  also  vary  from  cell  to  cell  affecting  histograms  and  moments  of 
data  gathered  from  a  number  of  range  and/or  azimuth  cells. 

The  examples  delineated  above  indicate  the  need  for  a  mathematical  model 
to  characterize  the  nonstationary  and  nonhomogeneity  of  radar  clutter.  Im¬ 
plied  in  the  examples  are  practical  radar  problems  such  as  cell  averaging 
CFAR  design  and  performance  evaluation,  accurate  prediction  of  first  order 
density  functions  of  clutter  power  with  attendant  calculation  of  false  alarm 
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and  detection  probabilities,  prediction  or  spatial  histograms,  etc.  In  order 
to  oe  useful,  the  model  should  have  a  firm  theoretical  basis,  be  relatively 
simple,  and  be  relacable  to  experimental  data. 


II.  APPROACH 


The  proposed  technique  for  characterizing  nonstationary  and/or  nonnomo- 
geneous  clutter  is  to  treat  the  first  order  probability  density  as  a  function 
of  random  parameters.  Thus  if  Z  is  the  instantaneous  value  of  clutter  power 
and  ap  &2>  • • •  are  parameters,  then  the  conditional  density  function  of  S 
is  d  ^  ( 3  Ax  ^ •••  °n) •  If  the  parameters  are  treated  as  random  variables, 
then  the  density  function  for  3  may  be  written  as  the  N  fold  integral: 

p(3)  •JJ*  ...  p1(2/a1>  a,,  ...  an)  p2(a1,  a2,  ...  aQ)  da1da2...dar  (1) 

where  p2(a^,  a?,  ...  a^)  is  the  joint  probability  density  of  the  random  pa¬ 
rameters  . 

Once  the  general  formulation  as  expressed  by  Equation  (1)  has  been  as¬ 
sumed,  the  functional  forms  for  and  P2  along  with  the  type  and  number  of 
parameters  (a^,  a2,  ...  aQ)  must  be  determined.  Several  possibilities  for  p^ 
exist;  however,  the  best  choice  appears  to  be  the  Nagakami  density  function, 
which  was  developed  to  model  the  fading  statistics  for  high  frequency  propa¬ 
gation. 

A.  Choice  of  Conditional  Density 

It  is  shown  in  [1]  that  that  Nagakami  density  is  an  approximate,  general 
solution  for  the  probability  density  function  of  the  quantity: 

r-  |L* 

i 

where  a^  and  ^  are  random  variables.  Note  that  Equation  (2)  can  be  viewed 


3^ 


(2) 
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as  che  instantaneous  envelope  of  clutter  which  is  modeled  as  the  sum  of  re¬ 
turns  from  indiviauai  discrete  scatterers  having  amplitudes  a.  and  phases 

L  a. 

Of  course,  che  instantaneous  clutter  power  s  is  just 


a  *  r 


The  functional  forms  of  the  Nagakami  density  functions  for  p(r)  and  p(z) 


p(r;k,u)  - 


,k  2k-l  -kr 
2(k)  r  a 

roc)  (u)k 


r  k  0 


p(»;k,u) 


,,  ,k  k-1  -ka 

<k>  «  aT 

rck)  (u)k 


■  k  o 


where  u  is  the  mean  or  average  value  of  a  •  r  ,  k  is  the  inverse  of  the  nor- 
ealizad  variance  of  a. 


Var  (a) 


and  T  denotes  che  gamma  function 
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f 


k  >  0 


It  is  seen  that  Equation  (4)  has  the  general  form  of  a  gamma  density  function 
where  the  parameter  k  is  a  positive  number  greater  than  zero.  For  the  Naka- 
gami  density  functions,  k  is  restricted  to  the  interval  k  >  -j*  either 
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ease  Equations  (3)  and  (4)  can  be  viewed  as  families  of  curves  with  different 
functional  forms  depending  on  the  parameter  k.  For  k  =*  1,  Equations  (3)  and 
(4)  become  Rayleigh  and  exponential  densities,  respectively.  For  k  = 

Equation  (3)  becomes  a  single  sided  Gaussian  density. 

As  derived  in  (1],  Equation  (4)  corresponds  to  the  Ricean  density  func¬ 
tion 


P(a) 


exp  ( — “) 

2a4 


exp( - r)  I 

202  0 


2a 


(6) 


2 

for  ao  ^  k  -  1.  In  Equation  (6)  a4  is  the  variance  of  the  random  quadrature 

2 

Gaussian  components,  A  is  the  steady  power  component  and  Iq  is  the  modified 
Bessel  function  of  the  first  kind  of  order  zero.  For  1  £  k  £  ,  Equation 

(4)  corresponds  to  the  Q  density  function. 


P(a) 


(—  +  -i) 

2  Sx 


v [!  (t 


(7) 


where 

a  *  2a2 (1  +  p) 

$  -  2a2 (1  -  p) 

2 

a  is  the  variance  of  the  quadrature  Gaussian  components  and  p  is  the  corre¬ 
lation  coefficient  between  the  quadrature  components. 

It  is  seen  from  the  above  discussion  that  the  Gamma  density  function  for 
the  instantaneous  clutter  power  as  given  by  Equation  (4)  is  justified  from 
theoretical  considerations  and  includes  the  exponential  and  Ricean  density 


functions  which  are  frequently  observed  in  experimental  clutter  data.  The 
Gamma  or  Makagami  density  function  also  simplifies  the  analysis  because  it  is 
a  function  of  only  two  parameters,  the  mean  (u)  and  the  inverse  of  the  nor¬ 
malized  variance  (k) .  Treating  the  two  parameters  as  random  variables,  the 
expression  for  the  instantaneous  power  becomes 

p(a)  *  II  p^Ca/u,  k)  p2(u,  k)  dudk  (8) 

_  jca 

.  ,,  ,k  k-1  u 

p(a)  m  f  -  p7(u,k)  dudk  (9) 

JJ  r(k)(u)K 

where  p2(u,k)  is  the  joint  probability  density  of  u  and  k  and  the  integrals 
are  taken  over  the  full  ranges  of  allowable  values  of  u  and  k. 

B.  Examples 

The  model  for  the  first  order  density  function  as  expressed  by  Equation 
(9)  is  directly  applicable  to  the  short  and  long  time  variation  in  the  scat¬ 
tering  from  a  given  cell.  For  example,  suppose  that  it  is  desired  to  model 
the  terrain  scattered  return  over  an  interval  of  wind  velocities.  It  has 
been  observed  that  for  certain  vegetation  the  functional  forms  of  p^(a/a^,a2) 
and  the  mean  of  a  vary  with  the  wind  velocity.  Thus  if  the  Nakagami  density 
is -a  good  representation  of  p^(a/ct^,  ct2)  >  Chen  knowledge  of  k  and  u  as  a 
function  of  wind  velocity  would  determine  p^(a/u,k)  in  Equation  (9);  the  fre¬ 
quency  of  occurrence  of  the  various  wind  velocities  would  determine  p2(u,k). 
Another  single  cell  terrain  example  is  the  variation  of  the  average  power 
level  of  returns  from  vegetation  depending  on  its  moisture  content,  density , 
height,  etc.  In  this  case  p1(s)  is  generally  exponential  (k  »  1)  and  does 
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not  change  its  functional  fora.  Thus  only  the  various  values  of  the  mean  (u) 
and  their  frequency  of  occurrence  for  various  conditions  would  be  required  to 
determine  p(a)  from  Equation  (9). 

Similar  examples  of  clutter  variation  on  a  single  cell  basis  occur  for 
backscatter  from  the  ocean  surface.  A  low  resolution  radar  will  generally 
observe  a  first  order  density  function  that  is  exponential  but  has  a  time 
varying  mean  depending  on  sea  state,  wave  height,  wind  velocity  and  direction, 
etc.  This  case  would  be  treated  in  the  same  manner  as  the  previously  dis¬ 
cussed  case.  High  resolution  radars  recording  returns  from  a  single  cell  of 
the  sea  have  observed  "spiky"  returns,  which  are  attributed  to  the  backscat¬ 
ter  from  individual  waves,  interleaved  with  returns  which  have  an  exponential 
first  order  density.  In  this  case  the  overall  density  for  p(a)  could  be  ap¬ 
proximated  by  the  discrete  version  of  Equation  (8) : 

p(«)  *  (1  -  P)  p1(a/k1,u1)  +  P  p2(a/k2,u2)  (10) 

where  P  is  the  probability  that  a  wave  is  in  the  cell,  p^Ca/k^u^)  is  the 
probability  density  of  the  clutter  power  when  no  wave  is  in  the  cell,  and 
P2(a/k2,U2)  is  the  probability  density  of  the  clutter  power  when  a  wave  is  in 
the  cell. 

The  same  techniques  can  be  applied  to  histograms  obtained  from  samples 
taken  from  different  range  and/or  azimuth  cells.  Assume  that  samples  from  N 
independent  cells  are  obtained  and  that  the  density  associated  with  the  nth 
cell  is  p(a/un,k^).  The  overall  histogram  can  then  be  written  as 


p(a) 


1 1  p<»/vV 


ai) 
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Ill .  INITIAL  RESULTS 


A.  Analysis 

In  order  to  obtain  some  meaningful  results  during  the  initial  phase  of 
the  program,  it  was  decided  to  investigate  the  effects  on  the  first  order 
density,  p(a) ,  caused  only  by  variations  in  the  mean.  Hence  the  parameter  k 
was  assumed  to  be  constant  in  Equation  (9) .  The  first  order  density  can  then 
be  written  as 


(k)k  ak“  1  e  u 

p(»)  -  I  - r -  p  (u)  du 

Jo  r<k)(u)lc  2 


where  p2(u)  is  the  probability  density  of  the  mean  or  average  value.  From 
Equation  (13)  p(a)  was  determined  for  various  representative  functions,  p2(u). 
Selected  plots  were  made  to  show  the  general  effects  of  p2(u)  on  p(a).  The 
P2(u)  densities  which  were  investigated  included  the  Gamma,  uniform,  log  uni¬ 
form,  and  log  normal. 

1.  Uniform 

A  frequently  used  analytical  density  function  is  the  uniform  densi- 


p  (u)  * - 

2  u2  **  ui 


u2  -  U  -  u] 


P  2 (u)  »  0 


u  <  u^, u  >  u2 
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Tnen  Equation  (14)  is  substituted  into  Equation  (13), 


?(a) 


u„  -  u. 


/ 


2  ,,,k  :<- 1 

(k)  a  e 


rca 

u 


du 


uL  ?(k)  (u)' 


(15) 


lc3 

Through  the  change  in  variable  t  *  —  ,  Equation  (15)  becomes 

ka 


p(a) 


k-2  -t  . 
t  e  dt 


(16) 


The  Incomplete  Gamma  function,  which  is  tabulated  in  several  sources,  is  de¬ 
fined  as 


y(a,x)  »  f  t*"*  e**C  dt 
'  o 


(17) 


for  a  >  0.  Thus  Equation  (16)  can  be  written  as 


p(o) 


(u2  -  ux)  r(k)  *•  ux 


YOc-1,  j»)  -  Y(k-1,  jp) 

u2  J 


(18) 


for  k  >  1.  For  the  special  case  of  k  ■  1,  Equation  (16)  becomes 


In  terms  of  the  exponential  integral : 


E,(x)  = 

i 


/cc  -t 

V 

T 


Equation  (19)  can  be  expressed  as 


»<■>  ■  [E1  {fj  -  h  (^)j 


for  k  *  1. 

In  order  to  determine  the  effects  on  p(e)  of  the  spread  in  the  mean  u, 

u2  i 

values  of  p(a)  were  calculated  from  Equation  (21)  for  dynamic  ranges  —  »  — 

U1  a 

of  1,11,  2,  and  10.  These  values  are  listed  in  Table  1  and  plotted  in  Figure 

1.  For  comparison  purposes,  it  was  decided  to  plot  the  quantity  u2p(e)  ver- 

sus  the  normalized  variable  —  .  Another,  possibly  more  useful,  plot  would 

u2 

show  u  p(-)  versus  —  where  u  is  the  mean  of  u  and  is  related  to  u^  and  u2 

u 

through  the  relationships: 


2  1  +  a 


1  1  +  a 


Figure  1  and  the  data  in  Table  I  indicate  that  for  small  spreads  in  the 

tl  T  3 

mean  (a  large,  —  *  —  small) ,  p(a)  approaches,  the  function  •  As  the 

U1  a  u2 

spread  in  the  mean  u  increases,  the  function  of  u2p(s)  becomes;  more  concen¬ 
trated  near  the  origin.  These  results  are  to  be  expected,  because  as  a  ap¬ 
proaches  1.0,  u^  approaches  u2  and  p2  approaches  a  delta  function  centered 
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l"w«»»i  „ 


".> Vv'fjif  •j.'-VW  *~  v 


2 

a  *  0.9 

a=»  0.5 

a=  0.1 

a  0 

U2 

U2p(e) 

U2p(a) 

U2p(e) 

U2p(a) 

0 

1.054 

1.386 

2.558 

.1 

.948 

1.201 

1.782 

1.823 

.2 

.853 

1.041 

1.304 

1.223 

.3 

.768 

.903 

0.992 

.905 

.4 

.690 

.784 

0.776 

.702 

.5 

.622 

.681 

0.621 

.560 

.6 

.560 

.592 

0..504 

.454 

.7 

.504 

.515 

0.415 

.374 

.8 

.450 

.449 

0.345 

.311 

.9 

.408 

.391 

0.289 

.260 

1.0 

.367 

.341 

•0.244 

.219 

1.1 

.330 

.297 

.206 

.186 

1.2 

.297 

.260 

.176 

.158 

1.3 

.267 

.227 

.1505 

.135 

1.4 

.241 

.199 

M29 

.116 

1.5 

.217 

.174 

0.111 

.100 

2.0 

.128 

.090 

0.054 

.0489 

Table  III-l.  U^pCa)  As  Function  of  -jj-  For  Exponential 
p^Ca/u)  and  taiform  p£  Cu) 
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at  u  «  uv  The  resulting  expression  for  ?(e)  is  then 


?(*)  *  f 

J  o 


X  U 

e 


:(u  -  u,)  du 


As  a  approaches  zero,  u.,  approaches  zero  and  Equation  19  becomes 


.  /•»  -t 

-L  /  - — 

u2  '  —  c 
U2 


—  sT(— : ) 
2  1  u2 


The  expression  in  Equation  (22)  approaches  infinity  at 


0  and  has  the 


largest  negative  slope  for  all  values  of  a  between  zero  and  one. 


2.  Gamma  v 

Since  the  Gamma  density  encompasses  a  relatively  wide  family  of 
functions,  it  is  instructive  to  investigate  the  behavior  of  p(a)  when  p2(u) 
is  Gamma  distributed.  In  this  case 


P?(u) 


k_  k--l 

(kj)  2  «  2  a  “ 

-  k2 
T(k,)  (u)  1 


where  u  is  the  mean  of  u  and  k^  is  the  inverse  of  the  normalized  variance  of 

u.  Using  Equation  (13)  the  expression  for  p(z)  becomes 

/•*  .k2u  kz  . 

k  k2  k-1  /  “T  T 

(k)*  (k„)  «  1  /  e  U 


-.k  ,  2  k-l  I 

(k)  (k2)  e  1 

r(k)  r(k2) 


k-k2  +  i 
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k2u 

By  the  change  of  variable  x  3  ~r~  >  the  integral  in  Equation  (24)  can  be  writ¬ 


ten  as 


k2“  .  k. 


O  k~k2 +  1 

u 


du 


kak. 


k-k2  +  1 


d'X 


(25) 


Using  the  identity 


/ 


co 


-(X  - 


V+l 


„V+1 


dX 


V  V 


My> 


o  X 

in  Equation  (25)  and  substituting  the  resultant  into  Equation  (24),  the  ex¬ 
pression  for  p(a)  becomes 


k+k2  k+k2-2 


rrr 


p(«) 


2(kk2)  (a)  _  (y— ) 

r(k)  r(k2)  (u)  2 


k2+k2  ^2- 


(26) 


For  the  special  case  of  k  *  1,  p^(z/u)  is  exponential, 


p(a) 


k2+!  k2-l 


(2)(k2)  (a) 


V1 

r(k2)  (5)  2 


V 


4k2a 


u 


(27) 


The  density  function  given  by  Equation  (27)  is  especially  interesting  because 
it  can  be  integrated  in  closed  form  to  find  the  false  alarm  probability: 
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?.(a*)  -  j*  ?U)  da  3  i  -  j 


o(a)  da 


Substituting  Equation  (27)  into  equation  (28) 


P.  ■  1  - 


2(k2) 


V1 


Vi 

<3>  2  \2-l 


r(k2)  (u)  2 


Let  t  *  \  - *  »  chen 


41c2«* 


’f  * 1  ’  vT  / 
r(k2)  (2)  *  J 


A  V  2_ 

/  VT 


v 


(t)  dt 


Making  use  of  the  identity 


tV  Kv.l  (t)dt  -  2V“L  F(v)  -  x"  KV(X) 


Equation  (30)  becomes 


Pf(a*) 


2(k2a*)  “  K,_.2 


4k9a* 


r(k2)  (u)  2 


Thus  for  the  special  case  of  an  exponential  density  for  p., (e/u)  and  a  Gamma 
density  with  parameters  u  and  k2  for  p2(u),  the  false  alarm  probability  for  a 
given  power  threshold  of  «*  is  given  by  Equation  (32) . 

Figure  2  is  a  plot  of  Equation  (27)  for  values  of  k2  equal  to  1.5  and  3. 
Table  2  lists  the  values  used  in  the  plot.  The  behavior  is  similar  to  that 
shown  in  Figure  1.  As  k2  becomes  large  (k2  «  3) ,  the  variance  of  u  about  its 
mean  u  becomes  small  and  p(z)  approaches  the  exponential  density.  As  k2  be¬ 
comes  smaller,  larger  spread  of  the  variable  u  about  its  mean,  u  p(z)  attains 

a  higher  maximum  at  —  *  0  and  falls  off  more-  rapidly  as  a  function  of  -=■  for 

u 

small  values  of  the  argument.  These  observations  are  verified  by  the  limit¬ 
ing  forms  of  k2>  As  k2  approaches  infinity,  P2(u)  approaches  6(u-u)  andp(z) 

J5.  i 

approaches  e  u  .  As  k2  approaches  the  limiting  value  of  in  the  Nakagami 

u 

density  function. 


3.2 


a 

u 

Figure  IIX-2.  Nornalized  Probability  Density  Function  of  Clutter  Power 
for  Exponential  -  Gama. 
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It  is  interesting  to  note  that  an  independent  analysis  of  histograms  obtained 


from  range  swept  radar  data  [2]  indicated  that  p(a)  had  the  form 


0(3) 


c  e 


-b  V7 


(21) 


Va 


which  is  the  same  fora  as  Equation  (33) . 

3.  Log  Uniform 

Radar  backs carter  data  is  usually  obtained  in  the  form  of  K  log  s 
instead  of  direct  measurement  of  the  quantity  s.  Thus  it  is  useful  to  inves¬ 
tigate  density  functions  for  p^(u)  which  are  distributed  in  terms  of  log  u. 
Consider  the  density  function 


P(y> 


7z  “  7i 


y2  -  y  -  7i 


(35) 


where  y  »  10  log^Q  u.  The  density  for  u  then  has  the  form 


P2(u) 


u2  >  u  >  ux 


u  in 


u. 


(36) 


where  u2  ■  (10) 10  and  u^  »  (10) 10 .  The  density  function  given  by  Equation 
(36)  is  the  log  uniform  density  function. 

When  the  log  uniform  density  function  for  p2(u)  is  substituted  into 
Equation  (13) 
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p(a) 


r 

w  ii 


(k)k  g1"-1  e 


kg 

u 


u. 


du 


TOO  (u)k+1  In  (— ) 


(37) 


For  the  case  k 


p(a) 


In 


/u„ 


u 


du 


(38) 


Equation  (38)  can  be  readily  integrated  and  it  is  found  that 


g 

u„ 


g 

u. 


p(«) 


-  e 


a  u  V 


(39) 


Equation  (39)  is  plotted  in  Figure  3  from  the  computed  values  shown  in 
Table  3.  The  parameter  a  is  defined  in  the  same  manner  as  that  used  for  uni¬ 
form  density:  u^  ■  au^.  It  is  seen  that  the  behavior  of  p(z)  for  log  uni¬ 
form  p2(u)  is  similar  to  that  for  uniform  and  Gamma  p2(u).  For  small  spreads 
of  p2(u) ,  a  *  1.0,  p(a)  approaches  the  exponential  density.  Large  spreads  in 
P2(u),  a  -  0,  causes  u2p(a)  to  attain  large  peaks  at  —  »  0  and  to  fall  off 

rapidly  as  —  increases  for  small  values  of  —  . 

u2  u2 
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2 

a  ■  0.9- 

a=*  0.5 

a=  0.1 

U2p(e) 

U2p(e) 

U,p(e) 

0 

1.055 

1.44 

3.91 

.1 

.949 

1.24 

2.33 

.2 

.854 

1.07 

1.48 

.3 

.768 

.923 

1.00 

.4 

.691 

.797 

.708 

.5 

.622 

.689 

.521 

.6 

.560 

.595 

.395 

.7 

.504 

.515 

.308 

.8 

.453 

.446 

.244 

.9 

.408 

.387 

.196 

1.0 

.367 

.335 

.160 

1.5 

.217 

.167 

.0646 

3.0 

.045 

.0228 

.0072 

Table  III-3.  U^pCa)  As  Function  of  —  For  Exponential  p^(a/u) 
and  Log  Uniform  P2(u). 
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4.  Log  Normal 

Another  commonly  used  density  is  the  log  normal  density,  whereas 
the  name  implies,  the  log  of  the  variable  is  normally  distributed.  Thus  for 
y  *  10  ios1Qu, 


(v-v)2 

“  2a2 

y 


p(y)  -  t=t 


(40) 


V2*  a 


where  y  and  are  the  mean  and  variance  of  y,  respectively.  The  density 
function  for  u  is  then 


P2<u> 


■lnu  -  lnu. 4 

-  (  2  * 
2or 


u  \j2tr'  a 


«  i  u  i  0 


(41) 


wt ire  lnu  and  a  are  the  mean  and  variance  of  lnu,  respectively.  From  Equa¬ 
tion  (13)  the  density  function  for  the  clutter  power  becomes 


p(e)  - 


f 


GOk 


&  .  (lnu  -  lnu)' 

T  2a2 

e  _ 


\JlTp  a  TOc)  u1 


lcrt 


du 


(42) 


which  reduces  to 
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£ 

u 


( lnu  -  Inu)  ‘ 
9  a2 


e  e 


VI7  cr  u2 


du 


(43) 


for  k  ■  1.  Thus  far  attempts  to  obtain  closed  form  expressions  for  the  in¬ 
tegrals  in  Equations  (42)  and  (43)  have  not  been  successful,  although  series 
representations  for  the  two  equations  have  been  derived. 

Nakagami  [1]  states  without  proof  that  the  density  function  for  t  *  10 
lag^a,  when  u  is  log  normal  and  p^(a/u)  is  a  Gamma  density  is  given  by  the 
expression 

Pt(t)  «  F(T,  k,  To>  S ( T,  k,  to>  (44) 


where 


F(x,  k,  tq) 


2  (k)k 

n  T(k) 


exp 


2k(x-x  )• 

Q 

m 


k  e 


2(T-To) 

m 


S(x,k,  xq) 


_ m 

o2  +  4  aQ2  kQ 


exp 


2  a  2  k2 (1-Q)2 

o 

m2  +  4  a  2  kQ 

o 


To  "  10  lo®10u 

xq  *  E(10  log^u)  *  mean  value  of  10  log^g u 
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Variance  of  10  iog^u 


a  »  20  10g  a  «  3.686 

io  Ujlio. 

transroraing  co  che  variable  2  -  (10)  *  a  10  ,  it  is  found  thac 


p(a) 


10  f  10lng\ 

alnlO  ?t  \  lnlQ  j 


1  1 


1  {V  .%S 

'  Um '  e  m  exp 

m 

.  2  Oka2  -S-)  . 

a 


3  r(k.)  yj  1  +  ka2  ~ 


(45) 


where 


0* 

um 

Inu 


»  variance  of  lnu 
>  median  value  of  u  «  e 
*  mean  value  of  lnu 


lnu 


For  che  special  case  of  k  ■  1, 


exp 


p(«) 


a2  (1  -  f-)2 

ax 

2(1  +•  a2  f-) 
_  m 


u  \/l  +  a2  — 
0  V 


(46) 


Equation  (46)  is  plotted  in  Figure  4  as  a  function  of  —  for  »  0.5  and  1.5. 


The  quantity  plotted  is  u  p(e) .  The  values  used  in  the  plots  are  given  in 
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u 

Figure  III-4.  Normalized  Probability  Density  Function  of  Clutter  Power 
for  Exponential  -  Log  Normal. 
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a. 

9 

a  =*0.5 

9 

a“  =*  1.0 

u 

up  (2) 

up(e) 

0 

1.65 

4.48 

.01 

1.61 

4.10 

.1 

1.31 

2.13 

.2 

1.06 

1.26 

.3 

.866 

.845 

.4 

.718 

.607 

.5 

.602 

.457 

.6 

.510 

.355 

.7 

.435 

.283 

.8 

.374 

.229 

.9 

.323 

.188 

1.0 

.281 

.156 

1.1 

.246 

.131 

1.2 

.215 

.110 

1.3 

.190 

.0934 

1.4 

.168 

.0797 

1.5 

.149 

.0682 

2.0 

.0853 

.0329 

3.0 

.0319 

.00876 

4.0 

.0133 

.00254 

5.0 

.00583 

.000773 

-6 

1Q.0 

.000141 

2.63x10 

Table  III-4. 

up  (a)  As  Function  of  3  For  Ex. 

orential  p^(a/u) 

u 

and  Log  Normal  P2(u)* 
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■  *  'iiiii.ii  'nhlrr  I  I""  "  rm- . --T— - r— — , . , 


~  (7~  =»  0.5).  For  larger  spreads  in  u 
chat  ac  che  point  a  *  0, 


?(°)  -  — 
a 


o 

(c~  large) ,  Equation  (46)  indicates 


(50) 


Thus  u  p(z)  at  :  ■  0  becomes  large  as  a“  becomes  .Urge.  This  trend  is  veri- 

2 

fied  by  the  plot  for  O  •  1.5.  It  is  seen  that  the  general  behavior  of  p(a) 
for  log  normal  u  is  similar  to  that  for  the  other  p?(u)  densities  discussed 
earlier. 

Nakagami  [1]  indicates  that  when  p2(u)  is  log  normal  and  c!~  is  large, 
p(a)  approaches  a  log  normal  density  when  p^(z/u)  is  a  Gamma  density.  For 
k  *  he  states  that  p(*>  is  log  normal  when  a  >  10  which  corresponds  to 

e  O 

a  >  2.3;  for  k  •  1,  p(e)  is  said  to  be  log  normal  for  a  >  7.0  or  a  >  1.61. 

o 

The  log  normal  density  has  been  proposed  for  terrain  at  low  grazing  angles 
[3]. 


B.  Data  Evaluation 

A  limited  amount  of  experimental  data  was  visually  inspected  to  obtain 
the  appropriate  form  of  p2(u).  It  appears  that  the  variation  in  the  average 
value  of  the  backscatter  from  trees  at  frequencies  of  9.5  GHz,  16.5  GHz,  35 
GHz,  and  95  GHz  may  be  approximated  by  either  log  uniform  or  truncated  log 
normal  densities.  Extensive  reduction  of  experimental  data  for  the  feature 
types  of  interest  at  the  desired  frequencies,  depression  angles,  time  of 
year,  etc.  would  be  required  to  properly  determine  che  correct  functional 
form  of  P2(u) ♦ 
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IV .  SUMMARY 


The  need  for  a  statistical  model  to  describe  nonstationary  and  nonnomo- 
geneous  radar  clutter  returns  is  discussed  in  Section  I.  Examples  are  given 
which  indicate  that  when  the  variability  and  uncertainty  of  clutter  are 
taken  into  account,  nonstationaricv  and  nonhomogeneity  of  the  returns  are 
the  rules  instead  of  the  exceptions. 

In  Section  II  the  general  and  specific,  approaches  to  the  problem  are 
outlined.  The  general  approach  is  to  model  the  first  order  density  function 
of  the  instantaneous  power  of  the  clutter  return  as  the  expected  value  of  a 
conditional  density  which  is  a  function  of  random  parameters.  The  expected 
value  is  taken  over  the  full  ranges  of  all  the  random  parameters.  If  n 
parameters  are  involved,  then  the  conditional  density  as  a  function  of  each 
of  the  parameters  must  be  known  along  with  the  joint  density  of  the  param¬ 
eters.  An  n  fold  integral  is  then  required  to  obtain  the  first  order  densi¬ 
ty  function  of  the  instantaneous  clutter  power. 

It  is  shown  that  the  Nakagami  or  Gamma,  density  is  a  good  choice  for  the 
conditional  density  function  for  several  reasons.  It  has  valid  theoretical 
and  practical  justifications  and  i3  a  function  of  only  two  parameters:  the 
mean  (u)  and  the  inverse  of  the  normalised  variance  (k) .  The  Nakagami  den¬ 
sity  includes  the  exponential  and  Ricean  densities  as  special  cases;  its 
particular  functional  form  depends  on  the  value  of  k. 

The  analysis  in  Section  III  considers  the  case  where  the  conditional 
density  has  the  Nakagami  form  and  only  the  mean  u  is  a  random  variable;  the 
parameter  k  is  assumed  to  be  a  constant.  Although  somewhat  restrictive, 
variation  only  in  the  mean  does  occur  in  practice  and  it  also  simplifies 
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Che  analysis.  Four  separate  density  functions  for  the  mean  are  considered: 
uniform,  Gamma,  log  uniform,  and  log  normal.  General  expressions  are  de¬ 
rived  and  selected  plots  are  made  for  the  first  order  density  of  the  clutter 
power  (a)  for  each  distribution  of  the  mean. 

When  the  mean  is  uniformly  distributed,  p(s)  can  be  expressed,  see 
Equation  (18) ,  as  the  difference  between  two  incomplete  Gamma  functions  for 
k  >  1.  For  k  «  1,  exponential  density  for  p(«/u)  ,  p(a)  is  given  by  Equation 
(21)  as  the  difference  between  two  exponential  integrals. 

For  a  mean  having  a  Gamma  density,  the  general  expression  for  p(»)  a3 
indicated  by  Equation  (26)  is  a  function  of  u  (the  average  value  of  u)  and 
k0  (the  inverse  of  the  normalized  variance  of  u) ;  it  also  contains  as  a  fac¬ 
tor  a  modified  Bessel  function  of  the  second  kind  of  order  k2  -  k.  When  k  *  1, 
p(»)  can  be  integrated  in  closed  form  to  obtain  the  false  alarm  probability, 
which  Is  a  function  of  kj,  u  and  the  threshold  level  a*. 

A  log  uniform  density  for  u  whan  k  •  1  results  in  an  expression  for 
p(a),  Equation  (39),  which  is  proportional  to  the  difference  between  two  ax- 
ponentlals  and  inversely  proportional  to  a.  The  exponentials  are  functions 
of  cho  upper  and  lower  limits  of  u. 

When  the  meen  hea  a  log  normal  danslty  the  Integral  expression  for  p(e) 
is  difficult  to  evaluate  in  closed  form.  Soma  results  which  are  stated  with¬ 
out  proof  in  (1]  were  used  to  obtain  an  expression  for  p(a),  Equation  (45), 

2 

which  is  a  complicated  function  of  k,  u  (the  median  value  of  u)  and  a 

m 

(the  variance  of  lnu) , 

For  the  case  where  the  conditional  density  p^(*/u)  is  exponential,  the 
geneml  behavior  of  p(a)  for  all  of  the  assumed  functional  forms  of  p,(u) 
give  similar  results.  Whan  the  spread  in  the  mean  is  small,  the  density  for 
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p(a)  is  approximately  equal  to  p^(a/u)  evaluated  at  u  =*  u.  As  the  spread  in 
u  becomes  large,  the  function  u  p(a)  becomes  highly  concentrated  at  a  =  0. 

Preliminary  visual  inspection  of  a  limited  amount  of  experimental  data 
on  the  backscatter  from  trees  seems  to  indicate  that  the  density  for  the 
mean  clutter  power,  p?(u),  can  be  approximated  by  log  uniform  or  truncated 
log  normal  densities.  In  order  to  obtain  reliable  estimates  of  p^Cu),  his¬ 
tograms  of  experimental  data  on  the  average  backscatter  for  the  terrains  of 
interest  at  the  desired  frequencies,  depression  angles,  polarizations,  etc. 


are  required. 


V.  FUTURE  WORK 


Future  efforts  on  the  characterization  of  nonstationary  clutter  should 
include  theoretical  and  experimental  investigations.  Several  analytical 
problems  are  suggested  based  on  the  effort  to  date.  First  of  all  the  re¬ 
sults  given  ir.  this  report  for  a  mean  which  is  log  normally  distributed  are 
based  on  an  equation  in  a  technical  paper  that  is  stated  without  proof. 

Since  clutter  with  log  normal  mean  has  been  measured  experimentally,  an  ad¬ 
ditional  effort  should  be  made  to  obtain  a  rigorous  derivation.  Another  ap¬ 
propriate  theoretical  topic  would  be  to  perform  an.  analysis  for  the  parame¬ 
ter  k  in  the  Nakagami  density  function  similar  to  that  reported  here  for 
the  parameter  u.  The  approach  would  be  to  assume  various  standard  densities 
for  p2(k)  and  integrate  the  expression  p1(a/k)p2(k)  to  obtain  the  probabil¬ 
ity  density  function  for  the  instantaneous  clutter  power  p(a).  Other  theo¬ 
retical  topics  would  be  suggested  from  the  analysis  of  experimental  data. 

For  example,  if  the  data  indicate  that  the  assumed  models  for  p^a/u)  and 
P2(u)  are  inadequate,  then  additional  derivations  similar  to  those  reported 
here  would  be  required. 

In  the  area  of  experimental  data  analysis  the  task  would  include  the 
generation  of  histograms  to  determine  for  various  terrains  the  conditional 
density  and  the  densities  for  the  mean  and  normalized  variance.  The  best 
fit  of  a. Nakagami  density  to  the  measured  conditional  density  would  be  made 
to  determine  the  adequacy  of  the  present  model  of  p^(a/u,k).  If  the  Naka- 
gami  density  is  found  to  be  valid,  then  temporal  or  spatial  histograms  of 
the  data  would  be  compared  to  p(s)  as  predicted  through  the  use  of  p^(a/u,k) 
and  p2(u,k)  to  determine  the  validity  of  the  overall  model. 
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